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Abstract — Among the priority issues that the modern society has to be solved, include also the energy and environmental 
issues. The notion of control (process control) has expanded in recent years, encompassing new areas such as automatic 
control of quality, the data processing with decisional purpose for one strategic leadership, ensuring uninterrupted of the 
system maintainability and thus, security and viability of the entire ensemble. In this context are part and the simplified 
simulation methods of the energetic installations from the power plants. The paper expose, in the first part, the importance of 
the steam generator in the operation of power plants, presents the energy processes complexity, emphasizing the importance 
of their management and automation to increase energy efficiency of each link in the chain. Then, taking as its starting point 
the real operational aspects, the mathematical modeling, the simulation and automatic control of steam generator, the paper 
proposes the development of a mathematical model in absolute units and the simulation the operation of a convective heat 
exchange surface of the steam generator in steady and dynamic regime. 
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I. Introduction 

The paper presents a study for the application about the steam generators in terms of modeling and simulation of their 
operation in steady and dynamic regime in power plants. Treat modeling itself of steam generator: it is based on physical 
model of each heat exchange surfaces, simplifying assumptions adopted, equations for conservation of mass, energy and 
momentum and the heat transfer equations. Also is presented the model in relative units imposed by the literature - 
highlighting the entire set of mathematical equations the processes occurring in the installation, the coefficients given in the 
literature, and the disadvantages involved in its use. Then, as an alternative, is proposed the steam generator model in 
absolute units: processing of the basic equations, are highlighted development its advantages - eliminating the recalculation 
sequence of parameters at each step, eliminating the linearization around the stabilized operating point, nonlinear 
characteristics are included in the model, determining the coefficients of thermal heat exchange between agents and the 
exchange surfaces in real-time, and the disadvantages involved in its use. The theoretical model proposed, in absolute units, 
implemented in Matlab-Simulink, leads to the simulations of the operating in dynamic regime for all components of steam 
generator. The model provides information about the influence of each parameter - input and output (pressure, temperature, 
and flow for water supply, combustion gases) on monitored parameters, namely pressure and temperature live steam. The 
results from the simulation are compared with values measured for parameters in real operation. The aim was to provide 
information allowing analysis of the model of each heat exchange surfaces on the one hand and on the other hand time for 
calculation to be within the real time operation. The steam generator model can be completed with the modules of the 
turbine-generator group with their installations and internal services and of the electrical part of power plant. All these 
models of the stationary and dynamic regimes can be implemented in computer systems and leads to the creation of a system 
for the management and control in electrical and thermal energy supply. 

II. The importance of steam generators in power plant operation 

Considering the great technical and economic importance for energy, the dynamics of the specific processes of the most 
important energy aggregates - steam generator and steam turbine - has been the subject of numerous studies and research, 
starting with the dynamic analysis of a steam generator with natural circulation and continuing with the analysis of its 
optimal and suboptimal control, or of the plant assembly. The vast majority of models is designed specifically for control or 
proposes new techniques to solve. Therefore, global models are conveniently reduced to allow for dynamic analysis of the 
plant assembly with significant, generally accepted simplifications. For this purpose, the theory of linear control, which 
describes the state of the boiler using linear differential equations with constant coefficients, is sufficient, [2], [4]. 
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Steam generators (boilers) are water vaporization systems that use the heat generated by the combustion of a fuel or an 
industrial process (heat recovery boilers). In boilers used in power plants, the steam from the vaporization process is 
overheated at a temperature higher than the saturation temperature of saturated steam to raise its enthalpy. It distinguishes 
two main circuits: water-steam and air-combustion gases. The steam-water circuit of the generator consists of pipe systems 
immersed in the flue gas ducts. From a functional point of view, the following heat transfer surfaces of this circuit are 
distinguished: Economizer - increases the temperature of the supply water to a value close to the saturation value, the heat 
transfer between the water and the flue gases is convective; Vaporizer - it ensures the passage of water from the liquid phase 
to the saturated steam, the heat transfer being carried out predominantly by radiation; The Primary Superheater - realizes the 
overheating of the saturated steam up to the temperature level desirable, heat transfer convective - radiation and Intermediate 
Heater - occurs in conventional steam boilers and provides an increase in steam temperature already released in the turbine 
high pressure body. The manner in which the heat transfer is effected on the one hand and the level of temperature required 
by the water-steam agent, on the other hand, requires the way in which these heat exchange surfaces are located inside the 
channels of combustion gases, [2], [4]. 

Inside the vaporizer tubes phase change takes place, the heat transfer coefficients have large values. In these circumstances it 
is necessary to place the heat exchange surfaces into the furnace where the outside of the tubes, the heat transfer coefficients 
are also large - achieving heat transfer by radiation is predominant. For the primary, intermediate superheater, and 
economizer case, the heat transfer coefficients can be obtained inside the pipes are slightly lower than for vaporizer, the three 
areas are located in the convective generator zone. High power boilers energy groups, 510 or 1035 t/h are Benson type 
boilers with n-shaped building or tower. The boiler primarily presentation intended purpose of the study, namely the 
modeling of operating modes. Each power unit of 315MW is equipped with two steam boilers of 510 t/h identical type 
Benson, forced to pass single variable evaporation point, built in the form of n. 

III. Convective heat exchange surface presentation 

In the dynamic modeling and control of steam generator are used mainly three types of one-dimensional models for the 
analysis of fuel combustion stability at low loads, the dynamic analysis of air-gas path and predetermining spatial distribution 
of water-vapor properties - models for studying the control pressure, temperature, level, dynamic load variation. The models 
are mathematically expressed by a set of equations from differential equations to partial differential equations; their nature is 
depending on the modeled process complexity and the purpose of their use. Modeling and simulation of heat exchangers is 
an essential problem of any simulation program and especially of steam generators, detailing the following general concepts, 
[11]: segmentation, disconnection, the biphasic mixture and thermal transfer. The location and complexity of the processes in 
steam boilers steam boilers had required that the process can be structured in technological modules. 

Thus the economizer, a surface of the convective heat exchange, is mounted in the end of the heating surface of the steam 
generator, resulting in a saving of fuel corresponding to the heat recovery. Economizer heating surface is usually made from 
steel pipe or, more commonly, cast iron pipe with wings in the form of double coil, transversely spaced apart, forming 
economizer package supported by the support tube and connected to cylindrical collectors - inlet and outlet respectively. The 
economizer is placed on top of the convective circulation before the outlet gases to the air super heaters. Water passes 
through the economizer flue gas counter. The physical model of the economizer is shown in Fig. 1. 



Fig.1. The physical model of the economizer 

where: F wl , F w2 are the mass flow of water input-output, kg/s; F gl , F g2 are the mass flows of flue gas input-output, kg/s; Qgp - 
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the heat transmitted from flue gas to pipes, W; Q pw - the heat transmitted from sheet to water, W; t wi , t w2 - the temperature of 
the water input-output, C; t g i, tg2 the temperature of the flue gas input-output, °C; t g - temperature of the flue gases, °C; t p - 
average temperature pipes, °C; p wi , p w2 - the economizer inlet, respectively outlet water pressure, bar; h wb h w2 , h gi , h g2 - the 
enthalpy of the water inlet-outlet or flue gas, kJ/kg°C; M p - the metal mass, kg; V - the total volume free, m 3 ; A- water flow 
area,m 2 ; f f - friction coefficient; di - inside diameter of the pipe, m. 

The physical processes describing the economizer operation are: single-phase flow through pipes; heat transfer from flue 
gases to pipes; heat accumulation in the pipe material; heat transfer from the economizer pipes to the heating water. 

The simplifying hypotheses adopted in the mathematical model are, [11]: the economizer is without boiling; the variables of 
the model satisfy the physical principles; the economizer is considered to be made of a single package; the model is with 
concentrated parameters and contains only time, not spatial derivations; gravitational and acceleration losses are not 
important; heat transfer is predominantly by convection. 

With the hypotheses presented, the mathematical model is described by the: 

The impulse conservation equation is reducing to: 

P W 2=P W X~P F ll (!) 


The heat transfer equation, predominantly by convection, from the gases to the economizer pipes: 

Q KP =k KP s(t g -t p ) ( 2 ) 

with k gp - the convective heat transfer coefficient from the gases to the economizer pipes, W/mK; S - the heat exchange 
surface, m 2 ; t g - the average flue gases temperature, °C; t p - the average temperature of the pipes, °C. 

The heat transfer equation, predominantly by convection, from the economizer pipes to the water: 

Qpw= k pw S (t p -t w ) ( 3 ) 


with k pw - the convective heat transfer coefficient from the economizer pipes to the water, W/mK; S - the heat exchange 
surface, m 2 ; t s - the average water temperature, °C. 

Convective heat transfer coefficients are determined using the criteria equations given in the literature: 




\ 0,65 


•Pr 0 , 33 -C, C, 


(4) 


, „ A (w-d 

k pa =0,02b- -< . 

d\ v ) 


Pr 0,43 • C, -C, 


(5) 


Relationships in which thermodynamic properties are calculated in real time using MATLAB functions for the fluid achieved 
value by the pressure and temperature. 

The mass conservation equation: 

For the water 


d 

dz 


(KpJ = F w \ - F w 2 


and for the gases: 


(6) 
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j- T (v,p,)=F„-F , 2 m 

The equation of heat accumulation in the economizer pipes: 

M c —it )=Q -Q (8) 

P P j W/ ^ 8P ^ pw 

The energy conservation equation, taking into account the mass conservation 
For the water: 


Vp—(h) = Q +F t h.-F 0 h 

1 wr'w \ w / xppw w 1 w 1 w 2 i 


w 2 


V p c ) = k s(t —t )+ F x J , — F ~c J 0 (9) 

w / w w ^ \ w / pw V p w) H’l vvr vvl w2 w2 w2 y ' 


and for the gases: 


V p —ih )- F ,h ,-F ,h ,-Q 

g^g g/ gl gl g2 gl ^g 


V p c -^—(t )=F,c J , — F 0 c J ~ — k .sit — t ) (10) 

grg g j^_\g/ g l gig 1 g 2 g 2 g 2 g* \ g p) v ' 


The system of equations to determine unknown operational sizes, in relative units, is the same as the superheater, the 
corresponding coefficients being given by the literature, [8], [9]. 

In absolute units, the unknown sizes are the expression, [7]: 

For the water pressure at the exit of the economizer: 


Pw2=Pw 1 + 


1 




W A W 


A 


V J 


• 0.25 


( 11 ) 


For the average water temperature: 


= j [Q pw + h wl F wi - h w2 F w2 ]— - d t 

* wPwP pw 


( 12 ) 


For the water temperature at the exit of the economizer: 


t 0 = 2t — t , 

w 2 w w 1 


( 13 ) 


For the average gases temperature: 

h = \K F s i- h g 2 F g 2-Q g ]-T7^—dT (14) 

v gpg c pg 

For the gases temperature at the exit of the economizer: 

t g 2=2t g -t gl (15) 


IV. Simulating of the economizer 

As stated in the presentation of mathematical model equations of heat transfer models used in relative units, data from the 
literature were replaced by equations of heat transmission by convection. For both the convective heat transfer coefficient of 
combustion gas in the pipes and water pipes were prepared from the calculation block diagrams, Fig. 2. 
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Fig.2. Diagram calculation of the convective heat exchange coefficient 

MATLAB functions were created for Prandtl number, friction coefficient, thermal conductivity, kinematics viscosity and 
specific volume of water being determined by the pressure and temperature of the fluid attained, Fig 3. 



di 


c!5 
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invl 


Product] 



ff 


Fig.3. Diagram calculation of the friction coefficient 

From the point of view of the simulation of the operation of this heat exchange surface was analyzed at unitary variations of 
the input sizes on the outlet water pressure and the required water flow. 



Economizer 


Fig.4. Block diagram - simulating of the economizer operation 
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Time offset: 0 

Fig.5. Water pressure variation at the exit of 

THE ECONOMIZER WITH THE ENTRANCE 
TEMPERATURE 



Fig.6. Water pressure variation at the exit of 

THE ECONOMIZER WITH THE ENTRANCE PRESSURE 



Time offset: 0 

Fig.7. Water flow variation at the exit of the economiser 


Analyzing the dynamic behavior of the steam generator economizer, step variations of the sizes - water temperature entering 
the economizer, economical input water pressure, combustion gas inlet temperature - it is found that the balance is achieved 
by damped oscillations. Output water pressure is greatly influenced by inlet pressure. 

V. Conclusion 

The model in absolute units is the basis for simulating the operation of the dynamic steam generator, [7]. He presents the 
major advantage of eliminating the sequence of recalculation of parameters at each step - in the models in relative units the 
whole set of coefficients is determined and entered into the model according to the parameters achieved in a previous 
reference regime; eliminates linearization around the stabilized operating point, nonlinear features being included in the 
model; it finds its own non-linear function points and updates parameter values; allows the determination of heat exchange 
coefficients between thermal agents and real-time exchange surfaces, their values obtained by computing within the limits 
recommended by the literature; removes the empirically determined coefficients, including the calculated exchange 
coefficients, from the heat exchanges exchanged between the thermal agents; it is characterized, with all its advantages 
presented above, as a slower model - relative to the one in relative units - due to the iterative calculation. The model in 
absolute units can be the basis for an intelligent system for simulation in real time operation of a steam generator. It is 
recommended to develop on other energy equipment. This supposes mathematical models in absolute units for the turbine 
and the electric generator, simulations of their static and dynamic regimes operation. 
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